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Magnetic frustration in the context of pseudo-dipolar ionic disor- 
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Abstract. - We consider an alternative to the usual spin glass paradigm for disordered magnetism, 
consisting of the previously unstudied combination of frustrated magnetic interactions and pseudo- 
dipolar disorder in spin positions. We argue that this model represents a general limiting case for 
real systems as well as a realistic model for certain binary fluorides and oxides. Furthermore, it is 
of great relevance to the highly topical subjects of the Coulomb phase and 'charge ice'. We derive 
an analytical solution for the ground state phase diagram of a model system constructed in this 
paradigm and identify magnetic phases that remain either disordered or partially ordered even at 
zero temperature. These phases are of a hitherto unobserved type, but may be broadly classified 
as either 'spin liquids' or 'semi-spin liquids' in contrast to the usual spin glass or semi-spin glass. 
Numerical simulations are used to show that the spin liquid phase exhibits no spin glass transition 
at finite temperature, despite the combination of frustration and disorder. By mapping onto a 
model of uncoupled loops of Ising spins, we show that the magnetic structure factor of this phase 
acts, in the limit T — >■ 0, as a sensitive probe of the positional disorder correlations. We suggest 
that this result can be generalized to more complex systems, including experimental realizations 
of canonical spin glass models. 
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Introduction. — Disordered phases of matter may be 
categorized in terms of the relative importance of frus- 
trated interparticlc interactions and quenched positional 
disorder ('frustration plus disorder'). For example, if the 
frustrated interactions are magnetic and the positional dis- 
order is uncorrelated then one arrives at a recipe for a spin 
glass, a concept which underpins much thinking about the 
nature of disordered states 01. However, it seems reason- 
able to ask if the combination of frustration with correlated 
positional disorder can lead to distinct behaviour that is 
not encompassed in the usual spin glass paradigm. 

Pseudo-dipolar positional disorder is characterised by a 
two-particle correlation function in reciprocal space, 5(q), 
that decays like a dipole-dipole interaction. It contrasts 
strongly with conventional disorder, in which (7(q) decays 
like a screened Coulomb interaction (an exponentially cut- 
off power law). Spontaneously generated pseudo-dipolar 
disorder occurs in ice rule ferroelectrics [2] and a mag- 
netic equivalent has been observed in spin ice 13]. How- 



ever, to our knowledge, there has been no comprehensive 
theoretical study of the effect of pseudo-dipolar positional 
disorder on magnetic properties (although Villain has pre- 
viously discussed the topic in the context of insulating spin 
glasses pi). In particular, the effect of varying the rela- 
tive strengths of competing exchange interactions within a 
system possessing pseudo-dipolar positional disorder has 
remained an unsolved problem. Such a study could be rel- 
evant, either ideally, or approximately, to many real sys- 
tems. For example, a combination of ice-rules positional 
disorder with magnetic frustration could be realised exper- 
imentally in certain inverse spinel ferrites p] or by the fam- 
ily of fluoride pyrochlores, exemplified by CsNiCrFg. In 
these systems two ionic species (e.g. Ni^+, Cr-^^) are dis- 
tributed over the pyrochlore lattice and Anderson showed 
that by minimizing the Coulomb interaction the distribu- 
tion of ions should obey the ice rules (although the perfec- 
tion of the ice rules in these systems has not been exper- 
imentally determined). However it is not clear that spin 
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glass states are a generic property of this experimental 
class [6l[7]. Hence it is relevant to ask the questions, what 
kind of magnetic states should one ideally expect, and 
what are their experimental signatures? Here we consider 
the case of ideal ice-rules disorder. Such a study is not only 
of interest to illuminate an alternative to the traditional 
spin glass paradigm but is also relevant to the concept of 
the Coulomb phase [O] , a general consequence of pseudo- 
dipolar correlations. Furthermore, there is currently much 
interest in 'charge ice' [10-12 , in which mobile electrons 



can adopt ice-rule configurations, leading to fractional ex- 
Our system provides a counterpoint to the 
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electron system, in that has quenched or static charge dis- 
order, which influences dynamical spin degrees of freedom. 

The Model. — The model we study consists of equal 
numbers of two species of classical Heisenberg spins s^ and 
Sf,, randomly distributed across the pyrochlore lattice (a 
cubic array of corner-linked tetrahedra) and subject to the 
'ice rule' constraint that there are two spins of each type 
on every tetrahedron. The spins are coupled by Heisen- 
berg exchange parameters that take one of three possible 
values, Jq, Jb or Jab depending on whether the neighbour- 
ing spins are both of type a, both of type b or one of each 
type respectively. In the case that all spin lengths and 
exchange parameters are identical, the model reduces to 
either a pyrochlore lattice ferromagnet or antiferromagnet. 
While the former orders conventionally, the latter remains 
in a cooperative paramagnetic state down to T = 13 
and gives rise to a distinctive pinch-point magnetic struc- 
ture factor indicative of pseudo-dipolar spin-spin correla- 
tions l3lfl4|H7l 



The case where spin and exchange parameters differ is 
addressed by expressing the spin Hamiltonian as a sum 
over contributions from individual tetrahedra. 
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(1) 



Here iV^ is the number of tetrahedra and {i,j)a indicates 
all pairs of spins i and j on plaquette (tetrahedron) a with 
exchange interaction J"-. 

A single tetrahedron. We begin by focusing on a single 
tetrahedron and, with no loss of generality, we assign spins 
Si and S2 to be of type a and spins S3 and S4 to be of type 
b. We then define three angular variables: (pa as the angle 
between spins Si and S2; (j)b as the angle between spins S3 
and S4; 9 as the angle between the resultants Sa = (S1+S2) 
and Sfc = (S3 -I- S4). In these coordinates the Hamiltonian 
for the single tetrahedron becomes 

H = -Jas1cOiicf)a- Jbsl cos (j)b 

-2JabSaSb{l +COS(/>a)'(l +COS0f,)5 COS 6*. (2) 

The ground state magnetic configurations then separate 
into two classes, characterized by the sign of Jab, allowing 
an effective phase diagram to be mapped out in the space 
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Fig. 1: The four regions of distinct magnetic behaviour 
(boundaries in white) as a function of the reduced exchange 
parameters Ja/\Jab\ and Jb/\Jab\- The boundary defining 
region IV is universal, |Ja||Ji,| = |JabP, while the remain- 
ing boundaries are defined by the ratios of the spin lengths: 
Ja/\Jab\ = —Sb/sa between regions I and II; Jb/\Jab\ = —Sa/si, 
between regions I and III. The background colour is an indica- 
tion of the value of the ground state sublattice order parameter 
as defined by equation [5] The bottom panels are the magnetic 
neutron scattering structure factors, S{Q), in the difi'erent re- 
gions as discussed in the text ((I-IV displayed left to right, 
respectively), shown in the [h, h, l\ reciprocal space plane from 
-2 to -1-2 on each axis). For region II, scattering from all spins 
has been included whereas in region III only scattering from 
components of the b spins perpendicular to the common a axis 
is shown (the intensity in this third panel is enhanced slightly 
relative to other regions in order to reveal the structure of the 
diffuse scattering). The plots of S{Q) have been obtained from 
Monte Carlo simulations at the points (Ja, Jb, Jab) = (1, 1, 1), 
(-2, 2, 1), (2, -2, 1) and (-2, -2, 1) for regions I, II, III and IV 
respectively (all at T/\Jab\ ~ 0.1). 



spanned by the reduced variables Jal\Jab\ and Jbl\Jab\ 
(Fig. [T]) . The topology of the resulting phase diagram is 
independent of the sign of Jab, as are the equations of the 
boundaries between regions. 

In region I, spins of the same species are aligned parallel. 
Spins of different species are either parallel or anti-parallel 
depending on the sign of Jab- This behaviour extends even 
into the quadrants of J-space where one or both intra- 
species interactions are antiferromagnetic. The extent of 
this domination by Jab is governed by the ratio Sa/sb as in- 
dicated by the equations of the boundaries between region 
I and regions II and III (Figure fl]). 

Regions II and HI are identical on interchanging the 
labels a and 6. In the ground state, spins of one species are 
perfectly parallel whilst the antiferromagnetic coupling of 
the other species is frustrated by the coupling Jab- These 
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spins cant away from the colhnear axis defined by the first 
species through an angle equal to half that between spins 
of type a;, 

2 P s2 
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where {x, y} = {a, b} {{b, a}) in region II (III) and | J^l > 
Sy/sx- The single species sublattice order parameter is 
defined in these regions as 
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where N^^ is the number of spins of type x and the sum 
is over all spins of type x (our definition anticipates ex- 
tending the theory to a macroscopic lattice). Combined 
with (pi) we see that the ground state order parameter is 
inversely proportional to J^: 
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This order parameter is defined only in terms of the x 
spins as the y spins are perfectly ordered with respect to 
each other (parallel or anti-parallel depending on the sign 
of Jab) throughout regions II and III. Within regions I and 

ix) 

IV we simply define TOgs =1,0 respectively. 

In region IV the intra-species antiferromagnetic interac- 
tions dominate, leading to configurations with 4>a — i'b = 
TT. Thus the third term of (pi) is zero and the Hamilto- 
nian is independent oi 0. In this region, spins of different 
species are effectively decoupled. For Jab > the confiu- 
ence of the phase boundaries is particularly interesting as 
at this point the ordered ferromagnet becomes degenerate 
with the antiferromagnetic spin liquid. 

Extension to the macroscopic lattice. Consider now a 
macroscopic pyrochlore lattice having two A and two B 
ions per tetrahedron. A walker, starting at some ion of 
type A and following a path only through sites populated 
by A ions (without retracing its steps) will always return 
to its starting point. Furthermore, the path traced out will 
have no branches but will form a continuous closed loop 
containing an even number of lattice sites. The whole lat- 
tice is tiled with such loops - every lattice site belongs to 
one (and only one) of these closed, even membered, loops 
of spins of a single species (referred to simply as 'loops' 
from now on). The statistics of such loops have recently 
been discussed in the context of the magnetic Coulomb 
phase 18 . In regions I, II and III, A type loops and B 



type loops interact with each other via Jab which has the 
effect of imposing long range order. By contrast, in re- 
gion IV there is no coupling between loops, although spins 
within a given loop are perfectly antiferromagnetically or- 
dered with respect to each other. We consider this soup of 
uncoupled, closed, antiferromagnetic loops to be a novel 
spin-liquid like phase; frozen interactions within loops ex- 
ist within a framework of two mobile degrees of freedom 



per loop. There are no energy barriers to facilitate global 
spin freezing and so, in the absence of free energy barriers 
or dynamical constraints, there can be no spin glass tran- 
sition in region IV. This is in contrast to the Heisenberg 
pyrochlore antiferromagnet (HPAFM) with weak random 
and uncorrelated bond disorder [l9p0] to which a number 
of points in region IV are closely related. 

In all regions of J-space, the ground states of the single 
tetrahedron are robust to stacking, with no extra frustra- 
tion incurred. The phase diagram in Fig. IT] should there- 
fore be equally valid for the macroscopic lattice. 

The decoupling of loops in region IV requires either 
(/)" = TT or 0^ = TT on every tetrahedron a, both of which 
are true in the ground state. The individual loops may 
then be viewed as independent one-dimensional Heisen- 
berg chains, which may be arbitrarily long in the thermo- 
dynamic limit. For T > 0, </>" and (f>^ may differ from their 
ground state values due to the excitation of low energy 
spin waves. The Hamiltonian then regains its dependence 
on the variable which is likely to dramatically slow the 
dynamics. Our numerical simulations suggest, however, 
that this slowing does not amount to truly broken ergodic- 
ity. Relaxing the ice rules constraint on the ion placement 
prevents this dynamical slowing down by providing extra 
unconstrained magnetic degrees of freedom in plaquettes 
with all ions of the same species or ions in a 1:3 ratio. It is 
readily shown that 1:3 and 3:1 tetrahedra place fewer con- 
straints on the magnetic degrees of freedom than do the 
2:2 tetrahedra. This result may be understood intuitively 
given that tetrahedra with all spins of the same species 
are the least magnetically constrained of all. 

Numerical Simulations. — Numerical evidence from 
Monte Carlo simulations is in agreement with our theo- 
retical predictions. Simulations were performed on lat- 
tices with L^ cubic unit cells {L — 4, 7, corresponding to 
1024 spins and 5488 spins respectively) for which a short 
loop algorithm l2l] was used to generate ion configurations 
obeying the ice rules constraint. In all simulations we ob- 
served loops on all scales, from the smallest possible (six 
membered rings) up to loops spanning the system. With 
one exception (discussed below) the energy scale was de- 
fined by |Jah| and a single spin flip Metropolis algorithm 
was employed, with spin updates confined to a small solid 
angle. 

To investigate spin freezing, we recorded the Edwards- 
Anderson order parameter 



Qea 



1 ^ 

-Y 



(6) 



in the region 0.01 < T/\Jab\ < 0-1 with L = 7. We chose 
Sa — 3/2, Sly = 1, (corresponding to the magnitudes of 
the magnetic ions in CsNiCrFg) and focused on the point 
{Ja/\Jab\,Jb/\Jab\) = (-1.1,-1.1), ensuring three differ- 
ent bond contributions to the energy. We simulated both 
the completely random and ice rules constrained models 
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Fig. 2: Variations in the order parameter, m^^\ along lines 
in J-space denoted {Ja/\Jab\,Jb/\Jab\)- (a) along (2, Ji,/| Jai,|); 
(b) along {Ja/\Jab\,'2); (c) along {Ja/\Jab\,-'2); (d) along 
(—2, Jb/\Jab\)- The solid lines represent the theoretical ground 
state magnetization as given by equationlSlin regions II and III 
(and defined as 1 in region I and in region IV). The symbols 
represent the results of Monte Carlo simulations conducted at 
T/\Jab\ = 0.01. Line symbols +, x, * and closed circles corre- 
spond to Jab > 0. Open symbols refer to Jab < 0. 



as described above. At each temperature the systems were 
annealed in five steps from T/\Jab\ — 1 with 10^ MCS/s 
(Monte Carlo Steps per spin) (10^ MCS/s were used at 
and below T/\Jab\ = 0.03) for equilibration at each step. 
Data was recorded over 10^ (10^) MCS/s and averaged 
over ten disorder configurations. In the absence of the ice 
rules constraint qea is essentially zero at all the temper- 
atures studied. Imposing the constraint leads to a signif- 
icant slowing of the dynamics, but with 10^ MCS/s qea 
remains below 0.05 even at T/\Jab\ = 0.01. As already 
noted, this behaviour is in strong contrast with that of 
HPAFM with weak random bond disorder 
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Figure^shows the variation of the order parameter m^'^' 
along the indicated lines in J-space. The agreement be- 
tween the theoretical and numerical results is striking and 
this level of accuracy has been achieved with relatively 
small scale simulations (L = A, Sa = 3/2, S{, = 1 and 
MCS/s=10®). These results validate the extension of our 
analytical solution to the macroscopic lattice. 

For the macroscopic lattice, differences in the ground 
state behaviour between the regions should manifest them- 
selves in the magnetic structure factor, examples of which 
are shown in the lower panels of Figure [l] For region 1, 
long range coUinear order produces sharp Bragg peaks. 
In regions II and III, sharp peaks arise from the compo- 
nents of all spins along the pseudo-collinear axis (shown 
for region II of Figure fl]), however scattering from just the 
perpendicular component of the canted spins (shown for 
region III of Figure fTl) reveals interesting diffuse scatter- 




Fig. 3: The magnetic structure factor S{Q) for a system tiled 
with independent chains of antiferromagnetically coupled Ising 
spins, each chain with a randomly selected easy axis. This 
data was obtained using L — 7 with averaging over 10 easy 
axis configurations for each of 50 disorder configurations. 



ing that suggests spin liquid like correlations transverse 
to the ordered component. We suggest that this might 
be called a 'semi-spin liquid', in analogy with a semi-spin 
glass (il. In region IV S{Q) has the characteristic struc- 
ture factor of an algebraic spin liquid, with pinch points in- 
dicative of the pseudo-dipolar correlations, although with- 
out the same clarity observed previously in studies of the 
HPAFM 
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The origins of such correlations in our 
model are not trivial. Unlike the pure HPAFM, the ground 
state in region IV has spins which interact only within a 
single loop. For such a ground state to exhibit pinch- 
point scattering would indicate that dipolar correlations 
emerge purely as a consequence of the geometric distribu- 
tion of loops, as governed by the ion configuration. The 
spin-spin correlation function is then nothing more than 
the probability that the two spins are on the same loop. 
To confirm this assertion we examined a toy model rep- 
resentative of the ideal ground state of the system: Jab 
was set to zero and perfect Neel order was enforced within 
each loop. We assigned a randomly selected easy axis to 
each loop and measured S{Q) for the resulting configu- 
ration. This process was repeated for a number of sets 
of randomly chosen easy axes and the resulting structure 
factors averaged. The resolution was improved by averag- 
ing again over a number of disorder configurations. The 
results (Figure p| clearly show signs of dipolar spin-spin 
correlations emerging from this system of magnetically in- 
dependent Id chains. To obtain the data in this figure we 
averaged over 50 disorder configurations with L = 7. S{Q) 
was averaged over 10 easy-axis configurations per lattice. 

A direct consequence of the above observations is that 
the magnetic scattering is acting as a probe of the struc- 
tural disorder. There is however an inherent limit on the 
resolution of the pinch-point scattering pattern which can 
result from a single realisation of the quenched disorder 
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on a finite lattice. Unlike the pure HPAFM, for which the 
bow-tie pattern is well resolved even for relatively small 
lattices, the ice rules constrained binary pyrochlore de- 
scribed here has the lengths and spatial arrangement of 
its loops predetermined by a particular ion configuration. 
The pure HPAFM however allows for a dynamic inter- 
change of spins between loops, in effect sampling a large 
number of loop configurations with a corresponding in- 
crease in resolution. In practical terms, this should have 
a similar effect to averaging the quenched disorder over 
many equivalent ion configurations. 

Conclusions. — In conclusion, for the ideal model 
considered we have demonstrated the suppression of spin 
glass behaviour and the emergence of novel spin liquid and 
semi-spin liquid phases. It will be of interest to re-examine 
the magnetic behaviour of the fluoride pyrochlores in the 
light of this result. In particular, certain fluoride py- 
rochlores have been reported to show spin glass transi- 
tions ml however there is a growing body of evidence that 
these compounds do not form traditional spin glasses be- 
low the supposed freezing temperature (tIIS 23 . We con- 



clude that if these are true spin glass transitions, they 
must be a consequence of disorder or interactions beyond 
those considered here. At a more general level we have 
illustrated a counter example to the idea that geometric 
frustration and positional disorder must combine to gen- 
erate a spin glass, although we have not ruled out the pos- 
sibility that the spin liquid states we have identified may 
be highly sensitive to further quenched disorder of a dif- 
ferent character. Finally, our results have shed light on an 
issue of rather general importance that is pertinent to the 
interpretation of neutron scattering patterns of disordered 
magnets. Thus, quenched atomic or ionic disorder is gen- 
erally characterised by an energy scale much higher than 
that of the magnetic interactions, so the magnetic struc- 
ture factor should generally be affected by the structural 
disorder correlations, but the question is to what degree, 
and what does this signify? We have identified a limiting 
case where the magnetic and structural disorder correla- 
tions are closely connected, and this connection is easily 
comprehended in terms of the model of 'unfrustrated' spin 
loops we describe. The opposite limiting case, which one 
would expect to apply to an ideal spin glass, is where the 
frustration is purely magnetic in origin and there is no de- 
pendence of magnetic on structural correlations. Real sys- 
tems are likely to lie between these two limiting cases, so 
our results may be of some general relevance to the inter- 
pretation of neutron scattering patterns, including those 



of canonical spin glasses 24 and weakly disordered spin 
glass systems such as 'SCGO' [25] and Y2M02O7 [26112^]. 
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